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Abstract. We probe the physical mechanism behind the known phenomenon of 
power synchronization of two diode lasers that are mutually coupled via their delayed 
optical fields. In a diode laser, the amplitude; and the phase of the optical field 
are coupled by the so-called linewidth enhancement factor, a. In this work, we 
explore the role of optical phases of the electric fields in amplitude (and hence power) 
synchronization through a in such mutually delay-coupled diode laser systems. Our 
mimerical results show that the synchronization of optical phases drives the powers 
of lasers to synchronized death regimes. We also find that as a varies for different 
diode lasers, the system goes through a sequence of in-phase amplitude-death states. 
Within the windows between successive amplitude-death regions, the cross-correlation 
between the field amplitudes exhibits a universal power-law behaviour with respect to 
a. 
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1. Introduction 

Synchronization in coupled oscillators has been the subject of much attention because 
of its fundamental importance in many areas of science and technology, such as in the 
dynamics of lasers [T], electronic circuits [2], biological systems [3], and its application 
especially in communication |H El EJ d El E] . Synchronization is known to be caused 
by the interaction between two nonlinear oscillators as a result of coupling, and, 
depending on the coupling strength and the time-delay in coupling, various features of 
synchronization emerge. The system of two coupled lasers is known to be an excellent, 
experimentally realizable example of coupled nonlinear oscillators, and synchronization 
in this context refers to the phenomenon in which the intensities of two lasers have a well- 
defined relation at all times. Different types of relations between laser intensities describe 
different types of synchronization [10], e.g., g'enerafced synchronization when intensities 
are functionally related to each other, complete synchronization when intensities of the 
lasers become identical (as a particular case of the generalized one), etc. There is another 
type called phase synchronization, when intensities of the lasers are uncorrelated but 
the phase difference of their oscillations remains bounded. This phase synchronization is 
usually defined as locking of the phases 9i^2, \pdi — q02\ < constant [H], where p and q are 
integers. The phase here indicates the phase of oscillations of the laser intensity, and not 
the optical phase of the electrical field. The phase synchronization of the oscillations of 
the laser intensities is studied by using the analytical signal concept proposed by Gabor 
[T^ [TTl [T3] . There exists a large body of work on synchronization and communication 
with chaotic laser systems - see Ref. [H] for a review. 

For the study of the synchronization phenomenon, mutually delay-coupled diode 
lasers are suitable candidates because of the compactness, low cost, and durability of 
diode lasers. Different aspects of the complex dynamics of mutually delay-coupled 
diode laser system have been probed - see Refs. [151 IE] and references therein. The 
system provides a simple and powerful tool to unveil the collective behaviour within a 
wide range of control parameter space, spanned by the coupling strength and the time- 
delay in coupling. Among the collective behaviour, amplitude death can occur in which 
two identical or non-identical coupled oscillators drive each other to a fixed point and 
stop the oscillations [17J. We have recently done experimental and theoretical studies 
of amplitude death and phase synchronization of powers in a mutually delay-coupled 
diode laser system [TSl HH]; however, to the best of our knowledge, the underlying 
physical mechanism of the power synchronization is not yet understood clearly. Given 
the importance of the system for fundamental studies as well as applications, we wish 
to extend our work to the understanding of the role, if any, of the optical phases of the 
laser fields in the process of power synchronization, for fixed values of coupling strengths 
and delays in the region of interest. 

This work can be put in the context of general research involving control of nonlinear 
dynamical systems using time-delayed feedback - see Ref. [2U] and references therein. 
In the systems studied in optics, passive feedback from mirrors or external cavities or 
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resonators has been used for this purpose. A notable mention is the work by Schikora 
et al. [21] for control of unstable steady-states in a semiconductor laser, which pointed 
out the role of the optical phase in the feedback control scheme. Our interest is in the 
understanding of the dynamical features of delayed coupling between two active devices 
of diode lasers, in the periodic regime of system parameters. In a diode laser, the active 
material has a highly asymmetric gain profile. This bears consequences to the refractive 
index (real part of the susceptibility), which can be related to the gain (imaginary part 
of the susceptibility) by the Kramers-Kronig relations [22]. The increase of the gain 
in diode lasers by increasing the carrier density leads to a decrease of the refractive 
index. The strength of the coupling between gain and refractive index is described by 
a parameter a, known as amplitude-phase coupling or linewidth enhancement factor 
[^ [21]. The a-factor influences several fundamental aspects of all semiconductor 
lasers, such as the linewidth, the chirp under current modulation, the mode stabihty, 
the occurrence of filamentation in broad-area devices. The dynamics of semiconductor 
lasers is greatly influenced by the a-factor, and its role has been specifically probed for 
dynamical effects such as instability enhancement in semiconductor lasers with delayed 
feedback [25] and injection-locking in semiconductor laser amplifiers 

In a coupled diode laser system, an amplitude fluctuation in one laser leads to a 
carrier density fluctuation, and through a, a phase fluctuation in the same laser. A 
change in the relative phase leads to an amplitude change in the second laser and an 
accompanying change in its carrier density. Thus a natural question arises: How does 
the phase relation between laser fields emerge in the complex dynamics? Does it have 
a bearing on the amplitude synchronization of the lasers? What is the effect of tuning 
the a-factor on the correlated amplitude dynamics? In this work, using the standard 
theoretical model for the system, we explore numerically the role of optical phases of the 
electric fields in power synchronization in a mutually delay-coupled diode laser system. 

The measured output powers of the lasers, P = \E\^ = oi course do not explicitly 
depend on the phases of the optical fields E, but there is indeed a connection between 
the field amplitudes A and the field phases through a. Our numerical results show 
that the synchronization of optical phases drives the powers of lasers to the synchronized 
death regime for a fixed a. Moreover, the cross-correlation measure between the 
field amplitudes shows that the system exhibits a sequence of in-phase amplitude- 
death regimes as a is varied for different diode lasers. In between successive in-phase 
amplitude-death regions, the cross-correlation coefficient shows a scaling behaviour with 
respect to a. It exhibits a power law, representing the transition from the in-phase 
amplitude-death to anti-phase periodic oscillations, and vice versa. 

In this paper, the theoretical model and numerical details to describe the two 
mutually delay-coupled diode lasers are presented in section 2. Section 3 presents the 
temporal dynamics of amplitudes and phases of the optical fields in such coupled lasers 
exhibiting synchronized death. The correlation measure is used in section 4 to explore 
the behaviour of synchronization over a large range of the amplitude-phase coupling a. 
This section also contains the scaling of the cross- correlation with respect to a within 
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the windows between successive in-phase amplitude death regions. The summary of the 
work is presented in section 5. 



2. Theoretical model 

A single diode laser subjected to optical feedback from an external cavity can be 
modelled by a set of fundamental delay-differential equations which are known as Lang- 
Kobayashi (LK) equations [2^. These equations describe the time evolutions of the 
complex electric field E{t) of a single longitudinal mode and the carrier density N{t) 
(with the threshold value subtracted out) averaged over the laser medium. In order to 
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Figure 1. Two diode lasers coupled in face-to-face configuration (schematic). 

analyse the behaviour of two diode lasers coupled in face-to-face configuration (as shown 
in figure [T]), the LK equations can be written in a standard normalized form [T8| [T9] : 
dEi 



dt 

dt 
dE2 

dt 
dN2 



T 



dt 



(1 + ta)Ni{t)E,{t) + r] e-'^'^E2{t 
Ji-N,-[2Ni + l]\Ei{t)\\ 
(1 + ia)N2{t)E2{t) + T] e-'^^^Ei{t 
J2-N2-[2N2 + l]\E2{t)\\ 



(1) 

(2) 
(3) 

(4) 



where rj is the coupling strength, i.e., the fraction of light of one laser injected into 
the other laser and vice versa, J's are the injected constant current densities (with the 
threshold value subtracted out), T is the ratio of the carrier lifetime to the photon 
lifetime, a is the linewidth enhancement factor as before, and r is the time taken by the 
light to cover the distance between the lasers. uji^2 are the optical angular frequencies 
of the solitary lasers 1 and 2. Ei^2{t — t) are the fields delayed by one coupling time 
T = L/c and Wi,2T are the phase mismatches. For simplicity, we have taken two identical 
diode lasers for which the detuning between them is assumed to be zero. In order to 
keep the model simple, we have ignored noise sources. The complex fields of the lasers 
can be written in terms of amplitude and phase parts as 

E, = A, e-*<^% (5) 

where Ai and 0, are the amplitude and the phase of the optical field of laser i {i = 1, 
2), respectively. Using ([5| in ([T])-(|4]), we get 

dAi 



dt 



Ni{t)Ai{t) + r]A2it - r) cos [0i(t) - - r) - Wsr] 



(6) 
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= J,-N,- [2N, + 1] \A^{t)f , (7) 
^ = - aN^it) - V^^j^ sin [Mt) - Mt - r) - cu^r] , (8) 

^ =iV2(t)A2(t) + r/Ai(t-r)cos[02(t)-0i(t-r)-u;ir], (9) 

dNo 9 

= J2 - - [2iV2 + 1] \A2{t)\' , (10) 

^ = - aN.it) - sin [Mt) -Mt-r)- cu.r] . (11) 

Note that the measured output powers of the lasers, Pi = {Eil"^ = Af do not exphcitly 
depend on the phases (pi of the optical fields. We still wish to probe the role of phases 
of optical fields in the power synchronization of two mutually delay-coupled diode lasers 
as a connects the phases and the amplitudes of the optical fields. Numerical integration 
of the above equations is done using Runge-Kutta fourth-order scheme with a step size 
= r/n, where n = 1000 is chosen based on the accuracy criteria. The dimensionless 
parameters are taken as Ji 2 = 0.165, and T = 1000 [IHl [IH]. It is found that the phase 
mismatch does not influence the results qualitatively and thus we keep ui^2T = — 1 (mod 
27r) [28j. Different kinds of lasers have different ranges of the linewidth enhancement 
factor a, e.g., for gas lasers, a is 0; for quantum dot lasers, a is 1.5 to 3 [291 EO], and 
for conventional diode lasers, a is typically 2 to 6 |3T] . 

3. Synchronized death 

We analyze the temporal behaviours of field amplitudes and optical phase difference 
in the amplitude-death regime, where the two mutually delay-coupled oscillators drive 
each other to a fixed point and stop the oscillations, as seen in [HI [19]. The temporal 
behaviour of amplitudes and phases of the optical fields are obtained numerically by 
integrating equations ([6])-(ll). In order to remove the initial transients, it is sufficient 



to discard the first 50000 data points before performing the synchronization study. We 
first choose a typical value of a = 5.6 corresponding to a diode laser, as in [TS] . 

From figure [2| it is clear that when the coupling rj between the lasers is zero, the 
laser field amplitudes (figure |2]^a)) and their optical phases (phase difference between 
the laser fields) (figure |2]^b)) do not synchronize to some constant values in a practical 
time scale. These are clearly shown in the corresponding insets, where the oscillations 
around 0.4062 (arbitrary unit) and —5.3292 rad, respectively, are still present even after 
a long time span. In this case, an amplitude change in one laser causes (through a) 
a phase change in the same laser, but the change in the relative phase does not lead 
to an amplitude change in the second laser. Thus due to the lack of this interaction, 
the amplitudes and the phases of the two laser fields do not synchronize to constant 
values. As the strength of coupling rj between lasers is switched on to a moderate value, 
e.g., rj = 0.2 and r = 14 (in units of cavity photon lifetime), the field amphtudes and 
optical phase difference synchronize to constant values of 0.7717 (arbitrary unit) and 
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Figure 2. The time (units of cavity plioton life-time) series of (a) the laser field amplitudes, Ai and 
A2 (continuous and dashed lines), and (b) the optical phase difference, </>i — </>2, for uncoupled lasers 
(r? = 0) at a = 5.6, with initial conditions Ai(0) 7^ A2(0), Ni{0) 7^ N2{0), and </>i(0) 7^ </i2(0). The 
initial transients are not shown. Insets show magnified temporal behaviours in the marked regions. 
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Figure 3. The same as figure [2] but for the coupled lasers with a coupling strength 7?=0.2, time 
delay r = 14 and a = 5.6. 



—3.1415 rad, as shown in figures |3](a) and |3[b), respectively. In the case when the 
couphng strength 7] between the lasers is non-zero, the optical phase difference of the 
fields approaches a constant value much earlier than the field amplitudes. This implies 
that the phase synchronization precedes amplitude-death synchronization, and it drives 
the field amplitudes to the synchronized death state. 

In order to quantify the decay of oscillations in laser field amplitudes and phase 
difference stated above, we analyze the variance ^32] as the measure of oscillations in 
the laser field amplitudes, Ai, A2 or the phase difference, 0i — 02, defined as 

a{Q) = AQ' = {{Q - {Q)f) , (12) 

where Q = A\^2 or 0i — 02, and (...) denotes time-averaging. This can be calculated by 
splitting the time series of field amplitudes and phase difference into non-overlapping 
time windows of duration At. After removing the transients, we calculate oiff) using a 
suitable At. We proceed to calculate the same for successive time windows (/ = 1, 2, . . .) 
of the same duration At, and stop when o < (0.001 x the value in the first time-window). 
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The values of ct{Q) are normalized to (t{Q), taking the maximum as = 1 at the 

first (/ = 1) time-window, and these data points are marked as circles in figures |4] and 
[5| for the uncoupled {rj = 0) and coupled {rj = 0.2, r = 14) cases, respectively. For the 
uncoupled case with a = 5.6, we choose At = 700 (units of cavity photon life-time). 
For the coupled lasers, the oscillations in field amplitudes and optical phase difference 
decay rather fast. Thus in order to generate sufficient number of data points, we split 
the time series of field amplitudes and optical phase difference into non-overlapping time 
windows of duration At of 350 and 112 (units of cavity photon life-time), respectively. 




2 4 6 8 10 12 2 4 6 8 10 12 



Time window index, / Time window index, / 



Figure 4. Scaled variance in (a) field amplitudes, Ai^2, and (b) optical phase difi'erence, (f>i — (p2, 
versus time-window index, Z, for the uncoupled system (r; = 0) at a = 5.6, using a time-window 
duration. At = 700 units. Circles are the numerical data, and continuous lines are fitted curves using 
equation ^ with (a) a = 2.296, -y At = 0.828 units, and (b) a = 2.330, 7 At = 0.844 units. 




Time window index, / Time window index, / 



Figure 5. Scaled variance versus time-window index, /, for (a) field amplitudes, ^1,2 1 using a 
time- window duration At = 350 units, and (b) optical phase difference, (/ii — 02i using a time- window 
duration. At = 112 units, for a coupled system with a time delay r = 14, a coupling strength rj = 0.2 
and an a = 5.6. Circles are the numerical data, and continuous lines are fitted curves using equation 
^ with (a) a = 3.562, 7 At = 1.270, and (b) a = 9.541, 7 At = 2.256. 
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The figures suggest that these oscillations decay exponentially, and can be fitted 
with the following expression: 

a{Q) = a e-^ ' ^\ (13) 

where a is a constant, 7 is the rate of decay per At, and / is the time- window index. 
From the fits in the uncoupled case {i] = 0), shown in figures |4]^a) and|4]^b), the decay 
rates (7 At) of oscillations in Ai^2 and (pi — 02 are found to be 0.828 and 0.844 per 
/, respectively. Thus, the oscillations in laser field amplitudes and in phase difference 
decay more or less with the same rate in the uncoupled case. 

Now, for the coupled system with coupling strength 77 = 0.2 and r = 14 units 



with a = 5.6, by fitting equation (13), we find the decay rates (7 At) of oscillations 
in Ai^2 and 0i — 02 to be 1.270 and 2.256 per /, respectively. Converting these decay 
rates per the same units of time-window duration as taken for the uncoupled case, 
we obtain (1.270 x 700/350) = 2.54 and (2.256 x 700/112) = 14.1 for the two cases, 
respectively, demonstrating that the oscillations in both field amplitudes and phase- 
difference decay much faster than those in the uncoupled case, and for the coupled 
system, the field-amplitude oscillations decay much slower than the oscillations in the 
optical phase difference between the laser fields. Thus, for the coupled system, the 
phase difference synchronizes fast, and it then leads to the synchronized death in field 
amplitudes. This is a central result of the paper, providing an insight into the mechanism 
of power synchronization in a mutually coupled diode laser system. 

4. Correlation coefficient 

The amplitude-phase coupling a describes the coupling between the real and imaginary 
parts of the susceptibility, and is given by the ratio of their derivatives with respect to 
the carrier density. Its value is known to depend on the carrier concentration, photon 
energy and operating temperature [33]. For gain-guided and low-dimensional lasers 
(quantum wells and quantum wires), the value of a can in fact be controlled by the 
design of the device structure. In order to analyse the role of a in the synchronization 
of mutually delay-coupled diode lasers of different kinds, we use a normalized cross- 
correlation function defined as 

^ ^ {{A,{t) - {A,m{A2{t) - {A^m) 

ViiMt) ~ {MtW){{A2{t) - {A^itwy 

We probe the characteristic features of C as the amplitude-phase coupling parameter 
a is tuned, and also mark the amplitude and phase dynamics at each distinct (a, C) 
point. The variation of C with a, in a typical range for diode lasers, is shown in figure 
[6| using the numerical solutions of equations ([6])-(ll ). The various fiat regions shown in 



figure |6] have a cross-correlation coefficient ~ 0.99. We have checked the dynamics and 
the cross-correlation values for higher values of the coupling rj and the delay time r. It 
has been seen that with higher values of 1], the fiat regions of amplitude-death shrink 
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Figure 6. Plot of cross-correlation C, given by equation | |l4| l, versus amplitude-phase coupling a 
for a time delay t = 14 (in units of cavity photon life time) and a coupling strength rj = 0.2. The 
different regimes are marked as X: in-phase amplitude death, and Y: anti-phase periodic oscillations. 



to smaller ranges of a, and with higher values of r, more transitions take place in the 
same range of a, the qualitative features remaining the same. 

In figure |6j the symbol X points to an in-phase amplitude-death region, while 
the symbol Y marks a region of anti-phase periodic oscillations. Corresponding to an 
in-phase amplitude-death state at a = 3.4 in figure [6} we show the synchronization 
dynamics in figure [7} For all such in-phase amplitude-death regions, field-phase 
synchronization precedes amplitude- death synchronization (as discussed in the previous 
section). This is again seen in figure [7| where the field-phase synchronization is 
established quite early, with the phase-difference 0i — 02 getting firmly locked at a 
constant value, as shown in figure [7]^c). Ai and A2 values are not exactly synchronized 
at these times, but the dynamics is shown to lead to the synchronized amplitude-death 
state at later times in figure [7|^b). The combined synchronization in the form of equal- 
time plot of the field amplitude-difference versus the field phase- difference, in the same 
time span, is presented in figure [7]^d). 

Within the window between two successive amplitude-death regions, say at symbol 
Y in figure |6] at a = 3.80, we have anti-phase periodic oscillations shown in figure 
|8j Here, the amplitudes and the phases of the laser fields are synchronized, as the 
oscillations in Ai — A^ and 0i — 02 remain within fixed bounds, around constant average 
values (figures |8](b) and (c)). The bounded synchronization is evidenced in figure |8]^d). 

In order to understand the transition from amplitude-death to periodic oscillations 
and vice versa, we explore one typical window between two successive flat in-phase 
amplitude-death regions in figure [6j The variation of the cross-correlation coefficient C 
with the amplitude-phase coupling a in such a window is shown in figure [9j Marked 
ttci and are the two critical values of a at which the cross-correlation coefficient C 
indicates transitions from in-phase amplitude-death to anti-phase periodic oscillations 
and vice versa. To characterize these transitions, we look for a scaling behaviour of C 
with a, as C varies sharply from 1 to —1 and —1 to 1 in figure |9} A power-law behaviour 
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Figure 7. Plots of (a) laser field amplitudes Ai and A2 versus time (in units of cavity photon 
lifetime), (b) magnified Ai — A2 versus the same time (with the inset showing further magnified 
behaviour in the marked region), (c) 0i — (p2 versus the same time (with the inset showing <j)i — (f>2 = 
constant even in a magnified scale in the marked region), and (d) equal-time Ai — A2 versus 4>i — <f>2 in 
the same time span, showing synchronization in the death state at a = 3.40, for a time delay t = 14 
units and a coupling strength r) = 0.2. 

is found as 

(1 + C) oc |ac - (15) 

where /x is the scahng exponent. The fitting of this scahng relation, shown in figures 
[IO|^a) and [lO|^b), gives the exponents jj, = 1.6 and 0.95, with ac ~ 3.7725 and 3.9630, 
respectively, for the two transitions. As mentioned before, a is different for different 
laser structures, and can be tuned by the carrier density or the light wavelength or 
the operating temperature to observe this critical behaviour. This general analysis of 
the complex dynamics can be extended to higher values of a, say, 4 to 11, suitable for 
quantum well devices 

5. Conclusions 

We have explored the role of optical phase dynamics in amplitude and power 
synchronization in a mutually delay-coupled diode laser system. It is found that the 
coupling a between the amplitude and the phase of laser fields indeed plays an important 
role in synchronization of such systems. The optical phase synchronization precedes 
amplitude-death synchronization, and it drives the field amplitudes to the synchronized 




Figure 8. Plots of (a) laser field amplitudes Ai and A2 (continuous and dashed lines) versus time 
(in units of cavity photon lifetime), (b) magnified Ai — A2 versus the same time, (c) ipi — 02 versus 
the same time, and (d) Ai — A2 versus rpi — <p2 in the same time span, showing synchronization in the 
anti-phase periodic state at o = 3.8, for a time delay t = 14 units and a coupling strength rj = 0.2. 
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Figure 9. C versus expanded a for one of the windows between successive amplitude-death regions 
shown in figure|6] «ci and a^re the critical values of a at which the transitions occur from in-phase 
amplitude-death to anti-phase periodic oscillations, and vice versa. 
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Figure 10. Plot of ln(l + C) versus ln|Q: — Qcj 2 I demonstrating the scaling behaviour of C with 
respect to (a) ctci and (b) ckc2 of figure |9] Dots are the numerical data and the solid lines are straight 
line fittings of the power law |15| . The critical values of the amplitude-phase coupling are found to 
be Qci ^ 3.7725 and a^a ~ 3.9630. 



death state. The approach to phase synchronization of the oscillating powers of the two 
coupled lasers can be easily observed experimentally at the appropriate time scale. The 
phase dynamics of a laser field, on the other hand, cannot be measured directly in an 
experiment - the phase difference between the fields can be revealed in the interference 
of the two fields. Our numerical prediction is expected to prompt experimental tests on 
such a system. 

We have also found that with the increase in amplitude-phase coupling a, the 
system hops over a sequence of in-phase amplitude-death regions. Within the windows 
between successive in-phase amplitude-death regions, the cross-correlation between the 
field amplitudes exhibits a power-law behaviour with respect to a, unveiling a remarkable 
universal feature in the dynamics of different kinds of diode lasers. 
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